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Abstract 

We investigate possible CP violating effects in e + e~ annihilation into top quark 
pairs. One of the interesting observable effects is the difference in production rates 
between the two CP conjugate polarized states tjj,^ and t^t^. The result is an asym- 
metry in the energy spectra of the lepton and the anti-lepton from the heavy quark 
decays. Another CP-odd observable is the up-down asymmetry of the leptons with 
respect to the reaction plane. These two asymmetries measure complementarily the 
absorptive and dispersive form factors of the electric dipole moment. Finally, as an 
illustration, we calculate the size of the CP violating form factors in a model where the 
CP nonconservation originates from the Yukawa couplings of a neutral Higgs boson. 
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I. Introduction. 

Since the top quark is widely believed to be within the reach of the present collider 
machines, it is not unreasonable for theorists to imagine what we can learn from the top 
quark. The best place to study the top quark in detail is in an e + e~ collider. One of the 
facts one would like to learn from the discovery of the top quark is the origin of the still 
mysterious CP violation. In this paper we investigate a way CP violation can manifest itself 
in the top pair production of an e + e~ collider. 

The top quark, due to its short lifetime, is believed to decay before it hadronizesjl]]. 
Therefore the information about its polarization may be preserved in its decay products. If 
that is the case, then one can investigate the source of CP nonconservation by measuring 
the CP violating observable involving a polarized top pair in the final state. This idea of 
detecting the rate asymmetry between different polarized states was recently proposed by 
Schmidt and Peskin[[|, |J. For ti production through the virtual photonic or Z intermediate 
states, to the lowest order in the final state quark mass, the polarizations of the quarks are 
either tjJ^B, or t^t^. (Note that we have adopted the notation that is the antiparticle of 

and should be left handed.) These two modes are CP self-conjugate. However since the 
top quark is heavy, there will also be large percentage of titi and £,r£r modes which are CP 
conjugates of each other. Therefore one can consider a CP asymmetry in the event rate , 
N(t L t L )-N(t R t R ). 

Schmidt and Peskin have shown how to detect the asymmetry N^iii) — N(tnif>) through 
the energy spectra of prompt leptons[§, ||. One assumes that the t quark decays semilep- 
tonically through the usual V — A weak interaction. Knowing that the hadronization time 
is much longer than the decay time|I|], one can analyze polarization dependence of its de- 
cay at the quark level. The top quark first decays into a b quark and a W + boson, which 
subsequently becomes l + v. For heavy top quark, the W + boson produced in top decay 
is predominantly longitudinal. Due to the V — A interaction, the b quark is preferentially 
produced with left-handed helicity. So the longitudinal W + boson is preferentially produced 
along the direction of the top quark polarization. Therefore the anti-lepton £ + produced in 
the W + decay is also preferentially in that direction. In the rest frame of the t, the angular 
distribution ||] of the produced £ + has the form 1 + cosip, with ip as the angle between £ + and 
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the helicity axis of the t. Above the ti threshold, the top quark is produced with nonzero 
momentum. As a result of the Lorentz boost, the anti-lepton £ + produced in the decay of 
the right handed top quark tn has a higher energy than that produced in the decay of the 
left handed top quark Similarly, the lepton I produced in the decay of II has a higher 
energy than that produced in the decay of t R . Consequently, in the decay of the pair t L t L 
the lepton from II has a higher energy than the anti-lepton from t^, while in the decay of 
tutu the anti-lepton has a higher energy Therefore one can observe N(tijti) — N(t^R) by 
measuring the energy asymmetry in the leptons. It turns out that this energy asymmetry is 
sensitive only to the absorptive parts of CP violating form factors. 

There is another equally interesting CP-odd effect in the azimuthal angular distribution, 
namely, the rate difference between the events with £ ± above the reaction plane and the 
events with below the reaction plane. Such an observable, like the previous one, is a 
direct measurement of CP violation and thus has no background from the CP conserving 
interactions. Unlike the previous case, this up-down asymmetry will probe the CP violating 
dispersive form factors. Though there have been many studies of CP violating observables in 
the literature p, |6|, we feel that the above two observables are simple and intuitive in nature 
and are easily implemented in future experiments. 

Finally, we will use a generic neutral Higgs model with CP violation in the Yukawa 
couplings to illustrate how these observables can arise. Such mechanisms of CP violation are 
contained in many extensions of Standard Model including the simple two doublet model [0. 
Throughout this paper, we focus our attention on CP non-conservation in the production 
mechanism only. There will be additional contributions if the usual V — A decay amplitude 
is also modified by CP violating interactions ||. 
II. CP violating form factors and amplitudes. 

We start by writing down the general form factors of the t quark. The vertex amplitude 
ieYi for the virtual 7* or Z* turning into t(p) and i(p') can be parametrized in the following 
expression: 

r £ = 4ln + 47m75 + + • • • , j = l,Z. (1) 
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We use the tree-level values for c v and c a . They are 

c 7 = - c 7 = 

c v = (i - fx W /)/y / x iy (l - Z W ) , (2) 
cf = -y^xw(l-xw) . 

Here xvi/ — 0.23 is the electroweak mixing angle in the Standard Model. The q terms are 
the electric dipole form factors. The spinor structure can be rewritten into another form 
using ijslp — p')^ = <J^(p + p'Y^s- Other irrelevant terms, like the magnetic moments, are 
not listed in Eq. ([!]). It can also be easily shown that, in the limit m e = 0, q is the only 
relevant form factor for the CP violating quantities we are interested in. 

The helicity amplitudes e 2 M(h e , h s , h t , hi) for the process e~e + — > ti at the scattering 
angle 9 have been given in the literature 0. For the initial configuration of cl&r, we have 

M(- + -+)= [c 7 + r L cf -/3r L cf](l + cos 6) 

M(- + +-)= [c 7 + r L cf + /3r L cf](l -cos^) 

M(- + — )= [2t(c 7 + r L cf)-|(c 7 + r L cf)/5/t]sin^ (3) 

M(- + ++) = [2t(c 7 + r L cf ) + |(c 7 + r L cf )/3/t] sinfl . 

Here we have used the convention B that CP invariance, when c^ z are turned off, is signified 
by the relation 

M(cr, a; A, A) = M(—a, — cr; —A, —A). (4) 

The dimensionless variables are defined by, t = m t /y/s, z = m z /y/s, (3 2 = 1 — At 2 . The 
Z-propagator and its coupling to the left-handed electron gives —er L /s with 

tl = {\- xw)/[0- ~ z 2 )\]x w (l - x w )] (5) 

The cross section is 

da(e L e R -> t ht t ht )/d{cos6) = §7ra 2 /3|M(-+, h t , h)\ 2 /s . (6) 

A color factor 3 has been included explicity in the above formula. Similarly, we obtain 
formulas for the initial configuration crCl with replaced by Tr, 

tr = -xw/[(l - z 2 )^Jx w (l - xw)] , (7) 
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and cos# by — cos#, and sin# by — sin# in Eq.(|3]). In case of an unpolarized e~e + machine, 
we must sum up initial configurations e^R and e R eL, and include the spin average factor ~. 

It is also interesting to note that, if the absorptive part of a scattering amplitude can 
be ignored (which is certainly true at the tree level), then the unitarity of the S matrix 
implies that the scattering matrix is hermitian. This hermiticity allows one to write down 
the constraint due to the CPT invariance as 

M(a, a; A, A) = M*(—a, -a; -A, -A). (8) 

We shall refered to this special case of CPT invariance as CPT invariance. The CPT invari- 
ance of course can be violated by the absorptive part of the scattering amplitude. 
III. Leptonic Energy asymmetry 

It is straightforward to see that the absorptive part Im or Im cf is required to produce 
the difference between configurations tjj^ and t R i R . This is expected because we need 
the absorptive part via the final state interactions to overcome the CPT constraint. The 
asymmetry integrated over the angular distribution is 



[N{t L t L ) ~ N(t R i R )] _ Z h=L , R 2P(cl + r h c z v ) (Im c] + rjm cf ) 
N(tt; all) " Z h=L J3 - + r h c z v Y + 2/3V£cf 2 " 



Note that there is no CP violating contribution due to cf coupling in the numerator when 
the electron mass m e is ignored. 

We can make use of this asymmetry parameter 5 to illustrate the the difference in the 
energy distributions of £ + or i~ from the t or i decays. The energy Eq{[ + ) distribution of a 
static t quark decay t — > £ + vb is very simpleQ in the narrow width approximation when 
rrib is negligible. It can be represented as 

{x (l — x )/D if mL/m? < x < 1, 

oi oj, w/t-o-, (1Q) 

otherwise. 

for on-shell W in the decay. Here we denote the scaling variable Xq = 2E (£ + ) / 'm t and the 
normalization factor D — g — |(mpi//m t ) 4 + \{mw I 'mt) 6 ■ When the t quark is not static, but 
moves at a speed f3 with helicity L or R, the distribution expression becomes a convolution, 

= /(so) 2xl(32 dx . (11) 



Here x = 2E(£ + )/E t . The kernel above is related to the (l±cos^>) distribution mentioned in 
the introduction. Similar distributions for the t decay are related by CP conjugation at the 
tree-level. Using the polarization asymmetry formula in Eq.(^), we can derive an expression 
for the difference in the energy distributions of t~ and l + : 

Ud£{Fj ~ 5^F)] = s \Jl(x,P) - h(xM . (12) 

Here distributions are compared at the same energy for the lepton and the anti-lepton, 
x(£~) = x(£ + ) = x = 4E(£ )/^/s. The count N includes events with prompt leptons or 
anti-leptons from the top pair production. It is useful to compare Eq.(|12|) with that of the 
overall energy distribution, 



1 r dN 



N 



+ 



dN 



dx(£+) dx(i-) 



Y, h=L ,RWr h c z a {cl + r h c%)[f R {x,P) - f L (x, 



+ h(x,P) + f R (x,l3) . (13) 

Here we only keep the dominant tree-level contribution. The first term of Eq. (|1^) is due to 

the two helicity modes that are CP self-conjugate. 

Fig. 1 
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Fig. 1 shows the overall prompt lepton energy distribution of Eq . (|13D , and the ratio of 
the expressions in Eq.(^) and Eq.(^) per unit Im (P d or Im cf . 



Note that our sample includes events with one or at least one prompt £ . However, if 
we are willing to use a smaller sample of prompt dilepton events (from the simultaneous 
semileptonic decays of both t and t) , a recent paper Ref . |10| has shown how one can use the 
kinematics of such events to analyze the helicities of t and i. Therefore such knowledge can 
be used to enhance the detectability of the asymmetry by cutting away the CP self-conjugate 
modes in the denominator of 5 in Eq.(|D |n|. However, the reduction in the size of event 



samples may be too high a price to pay. 
IV. CP— odd up— down asymmetry. 

It is known that explicit CP violation requires the CP nonconserving vertex as well as 
additional complex amplitudes. In the above case, this complex structure comes from the 
absorptive part due to the final state interactions. However, the complex structure can also 
come from other sources. One of these is the azimuthal phase exp(iL z (j)) in the decay process. 
This will produce a CP-odd up-down asymmetry even with only the dispersive part of CP 
violating vertex, Re q. 

The angular distribution of £ + from the t decay is specified by the the spin density matrix 
Pa,a' of the top quark. 

P {e)x,x' = M{ey 1 ]T M(h e ,h„x,ht)M*(h e ,he,x r ,ht) . (u) 

h e ,hi,hf 

Here M is the normalization such that Trp=l. p is hermitian by definition. 

dN(£ + ) = [(! + cos ip)p ++ + (1 - cos ^)p__ + 2Re (p+_e iv ) sin V>] dtpd cos V>/(4tt) . (15) 

The polar angles ip and the azimuthal angle <p of £ + are defined in the t rest frame F t , which 
is constructed by boosting the tt center of mass frame F along the motion of the top quark. 
In Fq, the z-axis is along the top momentum p t , the y axis is along p e - x p t , and the x-axis 
is given by the right-handed rule. The production plane is the x — z plane. 

Similarly, the angular distribution of l~ from the i decay is specified by the spin density 
matrix, Pxy, of the anti-top quark. 

p(0)x,y = E M(h e ,h S ,ht,\)M*(h e ,he,h t ,\') . (16) 

h e ,hg,ht 

dN(JT) = [(1 + cos^)p++ + (1 - cos^)p— - 2Re (p+_e~^) sin^]rf<prfcos^/(47r) . (17) 
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The polar angles ip and the azimuthal angle (p of i~ are defined in the t rest frame Ff, which 
is similarly constructed by boosting the ti center of mass frame F along the motion of the 
anti-top quark. It is important to keep in mind that the three coordinate systems Fq, F t 
and F[, have parallel directions of coordinate axes. 

In the phase convention^ such that Eqs.(|],[5[) are satisfied when CP is conserved, the 



following identities, first noticed by Gounaris et al.||ll||, 

p(0)x,x< = p(8)-x,-x< , (18) 

between the density matrix elements can be derived. Under CP conjugation, as we exchange 
i~ and £ + , by definition one should make the angular substitutions ip — > ir — ip, H> - > tt+V- In 
that case, one observes that the distribution in Eq.([T7|) is transformed into Eq.fll^D provided 
Eq.(|TBp is satisfied. On the other hand, the CPT invariance, when the absorptive amplitudes 
are ignored, implies 

p(o)x,x> = PW-k-x > (1 9 ) 

This is very similar to the CP and CPT transformations in the case of e~e + — > W~W + as 



analyzed before] 11, 12 



When the effect of the CP violating form factors are included in the analysis, one can 
form the following CP or CPT odd combinations: 

R(±)(B) X>X , = Re p(0) X)X , ±Re p(0) -X-X' ; 

and 

J(±)(0)a,a' = Im p(6) x ,x> ± Im p(6)-x,-x' ■ 

Among them, R(—) and /(+) are CPT odd; R(—) and 7(— ) are CP odd and all the others 
CP and CPT even. The observation of R(-) requires final state interactions due to CPT. 
Therefore it does not need to involve the the complex phase of the azimuthal dependence 



in Eqs.(15|, |17|). It can be decoded by analyzing the polar angular dependence of ij) or ip in 
these equations. In the collider C. M. frame these dependence can be translated into the 
energy dependence of the corresponding lepton in the final state, which has already been 
studied in the previous section. 
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Here we will focus on J(— ) which does not require final state interactions. Since p is 
hermitian, the only nonzero component of J(— ) is /(—)+_. It can be related to c d as 

N{9)lm \p(9)+- - p(6)- + ] = 2 sin^ £ (-/i)(c^ + r h cf ± (3r h c z cos^Re (c] + r h c z d )(3/t , 

(20) 

where the contributions h = R,L pick up the signs +, — respectively. One can in principle 
make detailed angular analysis of the difference between Eq . (|15D and its CP conjugate in 
Eq.(|17D similar to what was done for the case of e~e + — - > W~W + by Gounaris et al.fTT 



However, in an effort to find simpler observables which may be more intuitive and may 
be easier to detect, we shall consider the following partially integrated observable. Let 
dN(£ + , up) count events with £ + above the xz plane, i.e. p y {£ + ) > 0. Then, with other 
obvious notations, we define the following up-down asymmetry 

d = [dN(£+, up) + dN(£~, up)] - [dN(£+, down) + dNjT, down)] 
1 ' [dN(£+, up) + dN(£~,up)} + [dN(£+, down) + dN(£~, down)] ' 1 ' 

It is evaluated for each scattering angle 9. The branching fraction of the t semileptonic decay 
cancels in the ratio. Integrating on ip, (p or ip, (p over up or down hemispheres, we obtain 
A u ' d '(9) in a very simple form from Eqs. ([l5|JT?l) , 

A uA \9) = i(Im p(9)- + - Im P (9) + -) (22) 

As we sum up contributions from £ ± in each square bracket of Eq.fl2~ID, the asymmety is 
insensitive to the sign of charge, it is obvious that a non-vanishing value of A u d {9) is a 
genuine signal of CP violation. Although the angular distributions of the leptons derived 
from Eq.(^|) will have corrections from the strong interaction, the corrections cannot fake the 
CP asymmetry as the effects due to the strong interaction cancel away in the differences. 
To enhance statistics, it is useful to measure the integrated up-down asymmetry, 

A u.d. = I f ju.d.^ da d(cQS g) (23) 
a J a cosfcM 



From Eqs. (pO| - p3D , we obtain 



A uA. = _ Svr £fr=L,fl(-fr)(<ff + r h c z v )(Re cj + r fe Re c z d )(5/t 

16 E w (3-/5 2 )(c2 + r,cf)2 + 2/3V2 c f ' 1 } 
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Note that information on cf in Eq.fl20"D is lost when integrating the whole range of cos 8. 
However, one can easily find other convolution in Eq.fl2"3|) to recover the CP information due 



to ci 



Fig. 2 
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In Fig. 2, we plot the integrated up-down asymmetry w 4. n d - per unit Re or Re cf for 
various parameters. It seems that j\ u - d - per unit Re q increases with y^i, however, Re q is 
usually energy dependent, characterized by the underlying scale of new physics. Above that 
scale, Re q diminishes very fast. Therefore, the optimal choice of y/s is about at the scale 
of new physics. Note that even if Re Cd is relatively constant and the angular asymmetry 
is larger at higher energy, the event rate will become smaller because of the nature of the 
s-channel production. 

To measure the up-down asymmetry, we need a good determination of the reaction plane. 
This is possible for those events in which one of the top quarks decays hadronically into jets. 
One may be concerned about the imperfect angular resolution of hadronic jets. We argue 
that even if the angles of jets are ambiguous at the level of 10 degrees or so, as long as some 
sort of orientation of the reaction plane can be defined, the asymmetry will not be smeared 
away by more than one order of magnitude. For example, one can simply discard events 
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that the definition of up or down is made ambiguous by this smearing. The details of such 
a smearing effect will strongly depend on future detectors. 
V. Higgs Model. 

Among various mechanisms for the CP violation, the one that may manifest itself most 
easily is the neutral Higgs mediated CP violation. Since the neutral Higgs couplings are 
typically proportional to the quark mass, the large mass of the top quark naturally gives 
large couplings to the neutral Higgs bosons. CP non-conservation occurs in the complex 
Yukawa coupling, 

Ccpx = -(m t /v)i(AP L + A*P R )tH + {m 2 z /v)BHZ v Z v . (25) 

Here v = {\f2Gp)^^ — 246 GeV. The complex coefficient A is a combination of model- 
dependent mixing angles. Simultaneous presence of both the real part Ar = Re A and 
the imaginary part Aj = Im A guarantees CP asymmetry. For example, in the low energy 



regime, it can give rise to the electric dipole moment of elementary particles [JT3], [14]]. Here 
we will derive the CP violating form factors at high energy. They are induced at the one- 
loop level as shown in the Fig. 3. First, we calculate the absorptive parts according to the 
Cutkosky rules. The leading contribution to Imc^ comes from the rescattering of the top 
quark pair through the Higgs-boson exchange. Both CP violation and the final state effect 
are produced by the same one-loop graphs. 

Im ^O^(l4.o g(1 4)). ,26) 

The dimensionless variables are defined by, t = m t /^, z = mz/y/s, (3 2 = 1 — At 2 as before, 
as well as h = mnl J~s. For Imcf , there is a similar contribution. In addition, there could be 



a contribution due to the ZH intermediate state, Fig lb, provided the kinematics is allowed. 



cf T y _ aAjBc^t 2 
1 2(l-x w )x Vi 

Here /3| = 1 + h A + z 4 — 2z 2 — 2h 2 — 2h 2 z 2 , and the logarithmic factor 



Im cf = ^Im c] - " 1 y - [(3 Z + (2t 2 + 2t 2 h 2 - 2t 2 z 2 - h 2 )L] . (27) 



L-p^ l-f-V + Ph - (28 > 
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Our expression in Eq.(2?) agrees with that in Ref.[15]. Note that, at the threshold {(3 = 0) 



Im vanishes but Im cf has a value 



Im c 



z aA I Bc z u (3 z A-z 2 + h 2 (3 2 



8{l-x w )x w y -l~ z 2 -h 2 3(1 - z 2 - h 2 ) 2 ' ' 

In case mz + mu < \/s < 2m t , we still have the absorptive part Im cf , which is simply given 
by replacing the logarithm factor L in Eq. (^) by its continuation, 



2 Pz^TF 
arctan ■ 



(29) 



(30) 



l-z 2 -h 2 

The dispersive parts are obtained by the dispersion relation, 

u s 1 n f°° Im cj(s') , , 
Re cL(s) = -P / ^-^ds 

IX J s S f — S 

The symbol P denotes the principal value of the singular integral. Note that 

the absorptive part is evaluated at s' above the threshold sq, which is 4m 2 or 

(mz + ran) 1 for tt or ZH intermediate states respectively. It is understood that 

the dimensionless variables, (3, (3z, t 2 , h 2 and z 2 , are defined with respect to s'. 

Fig. 4 
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Fig. 4 shows typical sizes of d d for various cases. It is of order of 10 2 to 10 3 in gerneral. 
We must keep in mind that we only show the Higgs boson contribution in the perturbative 
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regime. In scenarios that the Higgs bosons interact strongly, the CP violating effect in the 
colliders could be much larger. 
VI. Conclusion. 

We have shown that the energy asymmetry or the up-down asymmetry of the prompt 

lepton events are sensitive to the CP violation in the top pair production in e + e~ annihilation. 

These two asymmetries are complementary to each other for measuring the absorptive and 

the dispersive parts of the CP violating form factors c^ z . 

Fig. 5 
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The cross section for e e + — ► tt at y/s = 400 GeV, for m t = 170 GeV, is ~ 1400 fb, 
among which about a quarter is due to the channels tjj,^ and t^pR (see Fig. 5). At an 



optimal luminosity of 10 cm -s , one year run of the futuristic NLCflni will produce 



20,000 prompt semi-leptonic events bhi^X. Note that we do not require both t and t to decay 
semileptonically. Therefore, only one power of the branching fraction 2/9 is involved. Such 
sample is good enough to analyze the CP violaing form factor a at a level of ten percent. 
Although predictions from Higgs models is still smaller than this level of sensitivity, the 
measurement could provide directly the important, model- in dependent information about 
the form factors of the top quark interactions. For a heavier top quark, y/s has to be raised 
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above the corresponding threshold, but not so far away as to incur the 1/s suppression in 
event rates. 
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Figures 

1. The energy distributions of prompt leptons, for the case that mt = 170 GeV, y/s = 
400 GeV. Case (a) for N~ l [dN / dx{i + ) + dN/dx(£~)], and case (b) for [dN/dx{i + ) - 
dN/dx(e-)]/[dN/dx(£ + ) + dN/dx(£-)} per unit of Imc] (solid) or Imcf (dashed). 

2. The weighted up-down asymmetry of prompt leptons, for the case that mt = 170 GeV, versus 
i/s, per unit of Rec^ (solid) or Reef (dashed). 

3. Feynman diagrams for the process e + e~ — > tt. The tree amplitude interferes with those one- 
loop amplitudes with (a) the final state interactions due to the exchange of a Higgs boson, 
or (b) the intermediate state of the ZH bosons. 

4. and cf versus yfs for the case that ran = 100 GeV and mt = 120 GeV. The parameters 
are chosen to be Aj = Ar = 5 = 1. 

5. The total cross sections are shown in the upper three curves for cases m^=120 GeV (solid), 
150 GeV (dashed), and 200 GeV (dash-dotted). The partial cross sections of the helicity 
configurations tjjj, and tRtR are shown in the lower three curves correspondingly. 
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Figure 3 
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